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PROJECT  REPORT  PA -229 -2 


The  channel  on  the  transcription  tape  labeled  Az  error  is  in  reality  traverse 
error.  Therefore,  please  replace  Paragraph  2  with  the  following: 

The  VHF  LC  RCS  is  computed  for  all  requested  range  gates.  The  gate  with 
the  peak  RCS  is  identified,  and  the  RCS  in  the  Az  and  El  error  channels  for  this  gate 
are  computed.  The  RCS  differences 


ATr  (db)  =  Az  error  channel  RCS  -  VHF  LC  RCS 
AE1  (db)  =  El  error  channel  RCS  -  VHF  LC  RCS 

are  computed  and  used  to  index  prestored  tables  of  ATr  (deg)  vs  ATr  (db)  and  AE1  (deg) 
vs  AE1  (db).  The  magnitude  of  the  angle  offsets  [ATr  (deg)  and  AE1  (deg)]  is  then  given. 
The  azimuth  offset  is  computed: 

AAz  (deg)  =  ATr  (deg)/cos  El  (deg). 


The  label  on  the  transcription  tape,  Az  error  channel  RCS,  has  not  been 
changed.  , 
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FOREWORD 


This  is  the  second  report  in  the  Data  Reduction  Program  Documentation 
series.  It  is  dated  according  to  the  date  of  completion  of  the  documentation.  No 
implication  is  made  that  this  program  will  not  subsequently  be  modified,  amended,  or 
superseded:  on  the  contrary,  the  history  of  radar  data  processing  is  one  of  continuous 
evolution  of  techniques,  and  it  is  unrealistic  to  assume  that  steady -state  has  been 
reached.  The  PA -229  series  is  being  published  for  the  convenience  of  interested 
parties,  and  Lincoln  assumes  no  responsibility  for  the  correctness  of  the  information 
presented,  nor  for  its  currency. 

The  preparation  of  reports  in  this  series  is  under  the  Editorship  of  Charles  R. 
Bemdtson  of  Lincoln,  and  of  D.  Nessman  and  R.  French  of  Philco-Ford  Corporation. 
Inquiries,  suggestions,  corrections,  criticisms,  and  requests  for  additional  copies 
should  directed  to  C.  R.  Bemdtson. 

The  principal  contributor  to  this  report  was  A.  J.  Poirier  (Philco-Ford).  Due 
to  the  intricate,  evolutionary  manner  in  which  the  programs  came  into  being,  the 
editors  regret  that  it  is  in  general  impossible  to  give  due  credit  to  all  --  mathematicians 
or  radar  analysts  or  programmers  --  who  contributed  to  the  definition  and  writing  of  the 
programs. 
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Alan  A.  Gromet^tein 
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ALTOAK 


PURPOSE  AND  UTILIZATION 


A.  Source  of  data 


ALTAIR1 


Data  Input 


v 

ALTAIR  transcription  tape 


w 

C.<  Description 

ALTOAK  is  used  to  analyze  data  onjargets  which  were  not  in  angle 


track.  For  any  object  in  range  track  by 


1,  the  program  computes  for 


selected  averaging  intervals  the  off-axis  angle  coordinates  (AAz,  AE1)  the  total 

ikjex*.  jmxAasv  tnt&b  .oecXwr-  s - a. .  ^ 

off-axis  angle  (A),  and  the^fteS-corrections  at  VHF.  and  UJdF.  The  angle  .  ^ 

coordinates  are  obtained  by  comparing  jthe  VHF^cT^S)with  the  RCS  in  the 

f  — i".  —"ij 

ancLfilerror  channels.  The  corresponding  phases  are  pfeed  to  determine  the 
s ens e  of  AA-g-  and^ABivy, _ 


D.  \  Outputs  y 

I  1.  A  listing  of  all  computed  quantities. 

'  2.  Plots  vs  TAL  of  the  uncdTrected  and  corrected  RCS  at  VHF 
LC  (Optional:  punched  cards  containing  thq/kCS  corrections). 

Punched  cards  containing  R,  Az,  and  El,  corrected  for 


known  errors,  in  a  format  Suitable  for, input  to  NRTPOD. 

\  /  / 


AtrT/vxR 

~#dA 


1 


CJJVC.V 


4  JUL-cJj-s^ 
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II. 


DESCRIPTION 


The  following  computations  are  performed  each  averaging  interval. 

The  VHF  LC  RCS  is  computed  for  all  requested  range  gates.  The  gate 
with  the  peak  RCS  is  identified,  and  the  RCS  in  the  Az  and  El  error  channels  for 
this  gate  are  computed.  The  RCS  differences 

AAz  (db)  =  Az  error  channel  RCS  -  VHF  LC  RCS 
AE1  (db)  =  El  error  channel  RCS  -  VHF  RC  RCS 

are  computed  and  used  to  index  prestored  tables  of  AAz  (deg)  vs  AAz  (db)  and  AE1 
(deg)  vs  AE1  (db).  The  magnitude  of  the  angle  offsets  (AAz  (deg)  and  AE1  (deg))  is 
then  given. 

The  pulse  by  pulse  phase  of  the  Az  and  El  error  channels  are  compared 
with  the  VHF  LC  phase  to  produce: 

N 

PHAZ  =  ^  ^  (Az  phase  -  VHF  LC  phase) 
i=l 

N 

PHEL  =  r,  /  (E1  Phase  -  VHF  LC  phase) 

N  Lj 

i=l 

where  N  is  the  number  of  pulses  in  the  averaging  interval. 

The  sign  of  AAz  (deg)  is  positive  when  cos  PHAZ  is  positive  and  the  sign  of 
AE1  (deg)  is  positive  when  cos  PHEL  is  positive. 

The  total  off-axis  angle  (e  )  is  found  by: 

9=  [  AEl^  +  AAz^  cos^  El  ]2 
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VHF  and  UHF  RCS  corrections  are  determined  from  tables  of  the  VHF  and 

UHF  beam  shapes.  Plots  of  corrected  and  uncorrected  VHF  LC  RCS  vs  TAL  are  pro 
duced. 


The  following  options  are  also  available: 

1.  Punched  cards  containing  TAL  and  the  VHF  and  UHF  RCS  cor¬ 
rections. 

2.  Punched  cards  in  a  format  acceptable  to  NRTPOD  may  be  obtained 
containing  R,  Az,  and  El  corrected  for  certain  errors.  R  is  corrected  for  bias, 
tropospheric  refraction,  and  target  position  in  sampling  pattern.  El  is  corrected  for 
bias,  tropospheric  refraction,  and  A  El.  Az  is  corrected  for  bias  and  AAz.  Ionospheric 
refraction  corrections  are  not  made. 

A  number  of  input  parameters  and  transcription  tape  parameters  are  checked 
for  validity  before  processing. 

The  main  program  checks  the  following  input  parameters: 

IP  AT  =  1  or  2 

TAVG  ±  0 

INTARG  f  0 

NRG  i  0 

2 

Subroutine  ALREAD  makes  a  number  of  other  checks  on  transcription  tape 
parameters.  For  some  errors  (missing  format  tables;  end  of  file;  target  no. , 
sampling  pattern,  or  polarization  not  on  tape)  information  is  returned  to  the  main  pro¬ 
gram  for  decision  to  terminate. 
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III. 


OPERATION 


A.  Input 

Start  and  stop  times  (GMT) 

Averaging  interval  and  skip  time* 

Target  and  sampling  pattern  numbers 
Specified  set  of  range  gates 

In  addition,  punched  card  output  may  b<"  requested  containing  RCS  cor¬ 
rections  or  observation  data  for  input  to  NRTPOD.  A  sample  ALTOAK  input  is  shown 
in  Appendix  A. 


CARD  1 

(15A4) 

(Col.) 

1-60 

TITLE 

60  character  title  for  printout  and  plots 

CARD  2 

(2(213,  F7.3),  4X,  415,  2F10.3,  215) 

(Col. ) 

1-3 

IH1  (13) 

| 

4-  6 

IM1  (13)  ^ 

>  Start  time  (GMT)  in  h,  min,  and  s 

7-13 

j 

ZSEC1  (F7.3)' 

14-16 

IH2  (13)  \ 

| 

17-19 

IM2  (13) 

i 

>  Stop  time  (GMT)  in  h,  min,  and  s 

20  -26 

ZSEC2  (F7.3)J 

1 

31  -35 

NRG 

Number  of  range  gates  (15) 

36-40 

INTARG 

Target  no.  (15) 

4 


41  -45 

IPAT 

Sampling  pattern  in  which  initial  gate  is  located  (15) 

46-50 

ING* ** 

Location  within  IPAT  of  initial  gate  (15) 

51  -60 

TAVG 

Averaging  interval  in  seconds  (F10.  3) 

61-70 

TSKIP 

Skip  time*  in  seconds  (F10.3) 

71  -75 

ICARD 

1  =  punch  cards  with  VHF/UHF  RCS  corrections 

0  =  no  punch 

76-80 

INRT 

1  =  punch  cards  for  NRTPOD 

0  =  no  punch 

CARD  3 

(315) 

If  INRT  =  0  Card  3  must  not  be  included 

1-  5 

IYEAR 

Last  two  digits  of  year 

6-10 

IMONTH 

Month  (1  to  12) 

11-15 

IDAY 

Day  of  month 

B. 

Output 

LISTING, 

GMT 

Az  and  El  corrected  for  bias 

VHF  LC  RCS  for  peak  gate 

Az  and  El  error  RCS  for  peak  gate 

AAz,  AE1,  and  9 

VHF  and  UHF  RCS  corrections  (db) 

PHAZ  and  PHEL 

CRLC  (corrected  VHF  LC  RCS) 

RFRANG,  AZCR,  RFELV  (corrected  values  of  R,  Az,  and 

El  for  input  to  NRTPOD) 


*Called  ISO  in  program  listing,  and  ISTGAT  in  ALREAD. 

**Skip  time  is  the  time  in  seconds  from  the  end  of  one  averaging  interval  to  the  start 
of  the  next. 


p;  ots 


Peak  gate  RCS  vs  TAL 

Peak  gate  RCS  corrected  for  off-axis  position  vs  TAL 
(Both  standardized  to  4  s/in  for  the  abscissa  and  20 
db/in  for  the  ordinate). 

PUNCHED  CARDS 


RCS  correction  data; 

TAL  (F10.3) 

AVHF  RCS  (db)  (F10.  3) 
AUHF  RCS  (db)  (F10. 3) 


Observation  data  (NRTPOD): 


Radar  identification  (A3) 

Year 

(3X,  12) 

Month 

(12) 

Day 

(121 

h 

(12) 

min 

(12) 

s 

(12) 

ms 

(15) 

Orbit 

no.  (IX,  I  1) 

Az 

(F8.3) 

El 

(4X,  F8.3) 

R 

(4X,  F12.4) 

R,  Az,  and  El  are  related  to  the  middle  pulse  of  an  averaging  interval.* 


They  are  determined  in  ALREAD  as  follows: 

R  =  R  +  R  (t  -  tR) 

0  0 

•  • 

where  R  ,  R  are  R  and  R  at  the  first  pulse  in  the  minor  cycle 

tp  is  the  time  of  the  first  pulse  in  the  minor  cycle 
t  is  the  time  of  the  middle  pulse 

Az  and  El  values  are  available  every  25  ms.  The  value  closest  to  t  in  the  major  cycle 
containing  the  pulse  is  used. 
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The  middle  pulse  is  determined  by  the  largest  integer  in  (N  +  l)/2,  where  N  is  the 

number  of  pulses  in  the  averaging  interval.  Other  output  quantities  are  associated  with 

the  median  time  of  the  averaging  interval,  determined  bv  (T  -  T  \/o 

last  pulse  first  pulse;/ 

bample  ALTOAK  outputs  are  shown  in  Appendix  B. 


IV. 


PROGRAM  LIMITATIONS 


A.  General  Remarks 

When  using  ALTOAK,  the  following  precautions  should  be  observed. 

The  noise  level  should  be  determined  theoretically  or  by  examining  a 
gate  that  does  not  contain  a  target.  The  signal  in  the  VHF  LC,  Az,  and  El  channels 
should  be  >  5  db  above  this  noise  level  to  obtain  valid  results.  If  the  error  channel 
signal  is  at  or  just  above  the  noise  level,  the  computed  angle  offset  is  an  upper  bound 
for  the  actual  offset. 

When  correcting  RCS  for  off-axis  position,  a  rule  of  thumb  frequently 
used  is  that  thn  cross  section  correction  must  be  s  6  db.  If  the  correction  is  >  6  db, 
the  accuracy  of  the  correction  is  questionable. 

The  user  should  be  careful  when  using  ALTOAK  for  trajectories  that 
have  high  Az  and  El  rates.  Due  to  the  method  used  in  ALREAD  for  picking  off  the  Az 
and  El  associated  with  a  pulse,  a  lag  of  as  much  as  50  ms  may  result.  This  lag,  which 
will  vary  in  a  saw  tooth  shape,  probably  will  be  small  compared  to  the  noise  in  the  Az 
and  El  offset  angles  for  trajectories  with  small  angular  rates.  The  RCS  corrections 
and  off-axis  position  data  are  not  affected. 


B.  Limits  of  Parameters 


Start  Time 

Stop  Time 

NRG 

TAVG 

TSKIP 

INTARG 

Length  of  Run 


Must  be  on  tape 
Must  be  on  tape 
30  gates 

Must  be  larger  than  the  PRI 
Can  not  be  negative 

Must  be  on  tape  within  start  and  stop  times 
^  1200  averaging  intervals 
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V. 


PROGRAMMING 


A.  OAKOS  (see  Appendices  C  and  D. ) 

OAK  OS  is  the  control  section  of  ALTOAK.  OAKOS  reads  the  input 
cards  and  calls  all  of  the  subroutines  that  process,  plot,  and  print  the  data. 

B.  GLMP  (see  Appendices  E  and  F. ) 

GLMP  searches  the  array  of  gates  selected  to  find  the  one  with  the 
laxg^st  RCS  for  each  time  interval.  GLMP  then  computes  and  saves  the  corresponding 
PHEL,  PHAZ,  Az  channel  RCS,  El  channel  RCS,  and  gate  number.  GLMP  then  inputs 
these  values  to  SUBOAK  along  with  GMT,  R,  Az,  and  lift-off  time. 

The  call  statement  is  GLMP  (AVGAL,  ISTGT,  AVGAZ,  AVGEL, 

IGAT,  AVGTM,  AVGRG,  NRG,  INDEX, 
ICARD,  TLIFT). 

INPUT 


AVGAL 

Alt* 

AVGAZ 

Az* 

AVGEL 

El* 

AVGTM 

Time*  (GMT) 

AVGRG 

R* 

ISTGT 

Start  gate  for  peak 

NRG 

Number  of  range  gates 

ICARD 

Flag  for  punched  cards 

TLIFT 

Lift-off  time 

OUTPUT 

INDEX 

Number  of  correction  c 

IGAT 

Peak  gate  number 

•Midpoint  in  averaging  interval. 


C.  TSPLIT  (see  Appendix  G. ) 


GLMP  calls  TSPLIT.  TSPLIT  is  used  to  convert  time  from  total  GMT 
to  h,  min,  s,  and  decimal  fractions  of  s. 

The  call  statement  is  TSPLIT  (AVGTM,  IHM,  TRUN). 

INPUT 

AVGTM  GMT  total  seconds 

OUTPUT 

IHM  (1)  Hours 

IHM  (2)  Minutes 

TRUN  Seconds  and  decimal  fractions  of  seconds 

D.  SUBOAK  (see  Appendices  H  and  J. ) 

GLMP  calls  SUBOAK.  SUBOAK  is  the  routine  in  which  most  of  the  com¬ 
putation  is  done.  It  consists  mainly  of  pre -stored  tables  and  equations  to  compute 
AAz,  AE1,  and  A  9.  SUBOAK  computes  correction  factors  for  UHF  and  VHF  RCS. 
STJBOAK  also  prints  the  data  and  punches  the  UHF/ VHF  correction  cards. 

The  call  statement  is  SUBOAK  (IHR,  IMIN,  ZSEC,  GMAX,  AVGAL, 

IGAT,  AVGTM,  AZE,  ELEX,  PHAZ, 
PHEL,  Az,  El,  INDEX,  AVGRG, 

ICARD,  TLIFT). 

INPUT 

IHR  Hours  (GMT) 

IMIN  Minutes  (GMT) 

ZSEC  Seconds  (GMT) 

GMAX  RCS  for  peak  gate 

AVGAL  Alt  for  midpoint  of  averaging  interval 

AVGTM  GMT  for  midpoint  of  averaging  interval 
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AVGRG 

IGAT 

AZE 

ELEX 

PHAZ 

PHEL 

Az 

El 

ICARD 

TUFT 


INDEX 


R  for  midpoint  of  averaging  interval 
Peak  gate  number 
Az 
El 

Az  phase  -  VHF  LC  phase 
El  phase  -  VHF  LC  phase 
Az  channel  RCS 
El  channel  RCS 

Flag  for  punching  correction  cards 
Lift-off  time 

OUTPUT 

Number  of  correction  cards  punched 


TIM 

TEM 

CRLC 

GTMAX 

RRANG 

RELEV 

AZCR 


STORED  IN  COMMON 
TAL 

Seconds  and  fractions  of  seconds  (portion  of  GMT) 

Corrected  VHF  LC  RCS 

Uncorrected  VHF  LC  RCS 

R  corrected  for  tropospheric  refraction 

El  corrected  for  tropospheric  refraction  and 
off-axis  position 

Az  corrected  for  axis  position 
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E.  REFC  (see  Appendix  K. ) 


SUBOAK  calls  REFC. 

The  tropospheric  refraction  correction  subroutine,  REFC,  is  based  on 

O 

tropospheric  refraction  tables  in  PPP-36.  A  modified  version  of  this  subroutine  is 
now  in  use. 

The  call  statement  is  REFC  (E,  R,  DEE,  DRR) 

E  Uncorrected  El  (must  be  between  0  ai  1  90  ) 

R  Uncorrected  R  (ft) 

DEE  El  tropospheric  correction 

DRR  R  tropospheric  correction  (ft) 

The  corrected  values  to  be  computed  after  exiting  from  the  REFC 

routine  are: 

El  =  E  -  DEE 

R  (ft)  =  R  -  DRR 

F.  OAKPLT 

OAKPLT  is  a  plot  routine  which  is  called  twice  from  OAKOS,  once  to 
plot  the  uncorrected  VHF  LC  RCS  vs  TAL  and  once  to  plot  the  corrected  VHF  LC  vs 
TAL.  The  user  has  no  control  over  the  size  of  the  plot. 


Q 


G.  REW 

4 

REW  is  an  entry  to  subroutine  BREADS  used  to  rewind  the  tape. 

H.  ALREAD2 

ALREAD  is  the  Fortrai  <  /er  for  the  assembler  language  tape  reading 
routines.  Appendix  L  presents  an  ALREAD  program  listing  which  has  minor  differences 
from  the  listing  in  Ref.  1.  These  differences  occur  because: 

1.  ALTOAK  requires  VHF  LC,  Az  and  El  error  channel  RCS 
data.  Only  one  channel  per  run  is  used  by  the  other  versions  of  ALREAD. 


2. 


Only  ALTOAK  requires  Az  and  El  data. 


The  call  statement  is  ALREAD  (TSTART,  TSTOP,  TLIFT,  INTARG, 
INPAT,  NOPHA,  NPTS,  NFP,  NEWPAS,  NRG,  JSTGAT). 


INPUT 


TSTART 

Start  time  of  processing  (GMT  total  seconds) 

TSTOP 

End  time  of  processing  (GMT  total  seconds) 

INTARG 

Target  number  to  be  processed 

INPAT* 

Sampling  pattern  in  which  initial  gate  is  located 

NOPHA 

2  (phase  and  RCS  data  wanted) 

NRG 

Number  of  range  gates  to  be  processed 

ISTGAT** 

Location  within  INPAT  of  initial  gate 

NPTS 


t 


INPUT  AND  OUTPUT  PARAMETERS 

Output:  number  of  pulses  of  data  returned 

Input:  must  be  initialized  by  calling  program 

before  each  call  to  ALREAD 


*Also  called  IP  AT. 

**Also  called  ING  or  ISG. 

+Set  to  zero  for  first  call.  Set  to  number  of  saved  points  for  subsequent  calls. 
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NEWPAS* 


TLIFT 

NFP 


Cycle  and  error  pointer  (see  error  returns  and 
switch  settings) 

OUTPUT 

Lift-off  time  (GMT  total  seconds) 

Frequency  code:  1  =  VHF;  2  =  UHF 


TIMES 

XSPH  A. 

RANGKM 

ALSAV** 

IRGA 

AZI 

ELE 

IPOL 

NPOL 


STORED  IN  COMMON 

Pulse  times  (CMT  total  seconds) 

RCS  and  phase  for  each  pulse  and  gate 
R 

Alt 

Range  gate  array  associated  with  XSPHA 
Az  (rad) 

El (rad) 

Data  channels  wanted 
Number  of  data  channels  used 


J.  Plotting  System  Subroutines 


The  subroutines  are  REREAD,  STOIDV,  and  PLTND. 


*Also  called  LAGAIN. 

** Valid  only  for  first  pulse  of  minor  cycle.  It  is  repeated  for  subsequent  pulses. 
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REFERENCES 


A LTAIR  Data  User's  Manual",  LM-97,  Lincoln  Laboratory,  M.I.T. 
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COMMON  SYMBOLS  AND  ABBREVIATIONS 


(The  units  given  for  certain  quantities  are  the  units  commonly  used  for  those 
quantities,  unless  otherwise  noted. ) 


ADT  ALCOR  Data  Tape 

Alt  Altitude  (km) 

APS  Average  Pulse  Shape 

ARS  ALTAIR  Recording  System 

Avg  Average,  Averaging 

Az  Azimuth  (deg? 

CADJ  Adjusted  Calibration  Constant  (db) 

C-band  ALCOR  frequency,  5664  MHz  (NB)  and  5667  MHz  (WB) 

El  Elevation  (deg) 

EOF  End  of  File 


GMT  Greenwich  Mean  Time 

h  Hours 

Hz  Hertz 


in  Inches 

LC  Left  Circular  Polarization 


t 


min 

Minutes 

NB 

Narrow  Band 

NRTPOD 

Non-real  Time  Precision  Orbit  Determination  Program 

POD 

Project  PRESS  Operation  and  Data  Summary  Report 

Phase 

Presented  in  deg 

PRF 

Pulse  Repetition  Frequency  (pps) 

PRI 

Pulse  Repetition  Interval  (s) 

pps 

Pulses  per  second 

pts 

Points 
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Range  (km) 

Range  Rate  (km/s) 

Radians 

Right  Circular  Polarization 
Radar  Cross  Section  (dbsm) 


Seconds 

Standard  Deviation  of  Wake  Velocity 
Time 

Time  After  Launch  (s) 

ALTAIR  Frequency;  415  MHz 
Velocity 

Doppler  Velocity 

Mean  Wake  Velocity 

ALTAIR  Frequency;  155.  5  MHz 

Wide  Band 

/ 

Tota .jl  Off-axis  Angle  (deg) 
Wavelength 

Denotes  Multiplication 
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FLOW  DIAGRAM  SYMBOLS 


PROCESS,  ANNOTATION 


DECISION 


> 


CD 


TERMINATOR 


SUBROUTINE:  where  NAME  is  the  entry 
call  into  the  subroutine 


CONNECTOR:  where  P  specifics  a  page  in  the 

flow  diagram,  and  L  designates 
a  statement  number  m  the  program 
listing  or  a  ’  -Terence  point  in  the 
flow  diagrf  n 

CONNECTOR:  where  X  implies  a  continuation 
of  the  diagram  to  the  next  page 


INPUT/OUTPUT  OPERATION 


MAGNETIC  TAPE 


c 


PUNCHED  CARD 


DISK 
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APPENDIX  B 
A  LTD  A  K  OUTPUT 


ALTA|»  DAK  VFR.MON  ?7  JIN  1  « 7 1 

ALTTAK  Cl  7  A I 
TARGET  NUMBFP  *  7 


T  1  *r 

1  44  3.?25 

l  44  3.2  75 

1  *4  1.32*5 

l  4  A  1.3  74 
l  <-4  i.4?S 

l  3.475 

1  44  3.5?5 

1  44  3.575 

1  44  3.624 

1  44  3.675 

1  44  1 .  7?  5 

1  44  3.775 

1  44  3.  F2  5 

1  44  3.9  74 

1  4  4  ->.Q25 

1  44  3,  <5  75 

1  44  4.025 

1  44  4.074 

1  44  4.124 

144  4.175 

1  44  4.225 

1  44  4.215 

1  44  4.326 

l  44  4.  3  75 

1  44  4.425 

1  44  4.4  75 

1  44  4.  525 

1  44  4.576 

1  44  4.626 

1  44  4.  475 

1  44  4.  724 

1  44  4.776 

1  44  4.P2S 

1  44  4  .  tj  75 

AT  T 1 HF  = 

1  44  4.C?5 


AT  T  IMF  » 


AT  T1MF  « 


n 

EL 

V-IC1C8 

50, 

10 

6, 

,re 

-  2 

•  58 

50, 

11 

6. 

►  07 

-2 

.21 

50. 

10 

6. 

,  f  7 

-2 

.40 

50, 

rq 

6. 

►  r  6 

-  1 

.87 

5  C . 

09 

6  , 

,06 

-2 

.28 

5  C« 

07 

6, 

f  5 

-2 

.  17 

50. 

05 

6. 

.04 

-2 

.48 

50. 

03 

6, 

.04 

-2 

.16 

50, 

00 

6, 

.03 

-2 

.17 

49. 

98 

6  , 

.0  2 

-1 

.86 

49. 

95 

6. 

02 

-1 

•  66 

49. 

92 

6 « 

.01 

-  1 

.25 

49. 

88 

6  , 

01 

-1 

.32 

45. 

94 

6. 

00 

-1 

11 

49. 

78 

5. 

99 

-0 

.90 

49. 

73 

5. 

99 

-0. 

,86 

49. 

68 

5. 

56 

-1 

.  29 

49. 

64 

5. 

,98 

-1 

.07 

49. 

53 

5. 

57 

-r , 

.  87 

49. 

47 

5. 

96 

-0, 

,  5  3 

49. 

40 

5. 

95 

-0 , 

.39 

45. 

’3 

5. 

55 

-0, 

.36 

49. 

26 

5. 

94 

-0 

.50 

49. 

22 

5. 

93 

-0  , 

.66 

49. 

14 

5. 

92 

-0, 

.*9 

49. 

07 

5. 

91 

- 1, 

,  18 

49. 

01 

5. 

90 

-0, 

.97 

49. 

57 

5. 

85 

- 1, 

.14 

48. 

92 

5. 

88 

-r , 

.84 

4  8. 

88 

5. 

87 

_  n , 

,9P 

48. 

P5 

5. 

F6 

-c, 

7  2 

48. 

80 

5  . 

P  5 

-0 , 

.53 

45. 

76 

5. 

84 

-0, 

87 

4  9. 

74 

5. 

84 

-1, 

.11 

6up.ci6c 

TARGET  7 

48. 

71 

5. 

P4 

-1. 

12 

fc?24. 

050f  THEPF  IS 

A  TIME 

GAP  08 

4 AO  DATA 

FKTEN 

fc22*  • 

1 0  Of  THERE  IS 

a  T  1ME 

GAP  OR 

BAD  CATA 

EXTFN 

Of  L7A 

DELTA 

DELTA 

DELTA 

DELTA 

AM  C.  I 

Fl  1081 

82  DIG 

H  DEC 

M6.EES  VWF(OB) 

UHFI08) 

-11. .4 

-14.01 

-0.70 

0.11 

0.89 

2.34 

13.37 

-n  .*6 

-4.81 

-0.61 

-C.  80 

1.09 

3.1  2 

27.87 

-IP. 71 

-12.27 

-0.76 

0.61 

0.99 

?  .8  7 

?  2  •  79 

-11.41 

-12.  17 

-0.62 

C.62 

0.88 

2.14 

15  .77 

-  1  1  .  B  1 

-  12.14 

-0.66 

C  •  65 

0.92 

?.  54 

17.20 

-11.11 

-4.  74 

-0.6  4 

-0.80 

1  .35 

3.16 

?8 , 53 

-12.41 

-11.76 

-0.«3 

C  •  1 6 

0.89 

?.  13 

14.47 

-12.42 

-1  1 .13 

-P.14 

-0.64 

0.90 

?.  44 

16.  46 

- 1  1.43 

-4.  72 

-0.61 

C.BO 

1.33 

3.04 

26.31 

-12.21 

-4.  74 

-0.14 

C.  77 

0.  97 

?.  78 

?  3.  2  8 

-12.04 

-13.42 

-0.14 

0.53 

0.79 

1.83 

1  7.93 

-  1 2 . 6P 

-4.  13 

-0.  IP 

-0.74 

0.91 

?  •  55 

17. ?7 

-13.14 

-  12.64 

-0.46 

C.  16 

r  .72 

1.35 

1  9.40 

-14.41 

-13.91 

-0.31 

0  .44 

9.6) 

0  .6? 

6.78 

-13.13 

-4.  66 

-0.  34 

-C.  72 

9.81 

1.07 

1  3  .67 

-13.41 

-4.74 

-0.36 

-0  .7  1 

0.89 

1.85 

13  01 

-11.44 

-4.  *  7 

-0.33 

-0.75 

0.8? 

2  .00 

13.78 

•11.61 

-4.14 

-0.32 

-C.  76 

0.  8? 

?.n4 

13. 9<J 

-14.76 

-13.63 

-0.27 

0.44 

9.  55 

0.3P 

6.06 

-1 1.41 

-4.  40 

-0.24 

-0.71 

9.77 

1  .67 

12.09 

-14.41 

-1.41 

•o.?i 

-f.  77 

0.  83 

2.  I' 

14.31 

-12.37 

-13.02 

-0.44 

-0.44 

0.65 

0  .96 

8.?? 

-1*.12 

-  10. fP 

-0.  3  4 

-C.63 

9.71 

1.31 

19.18 

-14.06 

-IP. 87 

0.11 

—  C  •  6  3 

9.7? 

1.  36 

1  9.  44 

-14.07 

-12.03 

0.31 

-0.15 

9.65 

O.Q? 

7.99 

-12.  16 

-  12.  44 

0.14 

-0.  56 

0.77 

1.71 

12.39 

-10.40 

-12.37 

0  .63 

-0.51 

0.83 

?.  r 

1  4.  31 

-1  2.  66 

-11.18 

0.49 

-6.64 

0.89 

1  .89 

1  3.23 

-11.24 

-10.74 

0.  19 

-C.65 

0.87 

2.  31 

15.M 

-1  1  .71 

-12.09 

P  .  16 

-0.17 

0.89 

1.87 

13.19 

-11.67 

-  10. 41 

0.  14 

-r.  63 

0.83 

?  ."7 

14.18 

-4,71 

-4.20 

0.69 

-C.  79 

I.05 

3.13 

?  8 .  r  1 

-11.11 

-14.14 

0.17 

0.42 

'-.71 

1  .78 

1 0 . 

-  10.  70 

- 10. r a 

0.  61 

-C.71 

0.96 

2.7? 

1  8.49 

IS  NOT  Pf<  THF  TAPE,  7 64 GETS  AVAILABLE 

ARF  1? 

-  1  n. 6  3 

-10.51 

0.66 

-r  .68 

0,75 

?  .  64 

17.92 

RFELV 

42CR 

f  ,*noo 


r. aoo3 


<>20.276 
OFO. 135 
979.863 
977.555 
979.250 
978.945 
979.638 
979.296 
977.989 

977.6*3 

°77.  375 
e77.033 
976.  727 
976. 42F 
976.090 
°75. 773 
c75.  466 

975.1?4 
9  74.  FI  A 
974.610 
974 .1 69 
973.  f  61 
973.  *64 
973.213 
972.906 
°72.600 
97? .294 
971.953 


9.227 
7.291 
7.  433 
7.429 
8.610 
7.330 
7  .466 
7.  T-*7 
7.  397 
3.452 
7.551 
7.  «45 
7.572 
7.  539 
7.48? 
7.402 
7.546 
7.5f  5 
7.454 
7.439 
7. 49  A 
7.490 
7.  4P8 
7.586 
7.467 
7.  556 
7.543 
7.410 


50.519 
50.431 
50. 548 
5^.546 
50,525 
50.477 

50.512 
50.451 
5". 454  . 
50 .477 
53.467 
50.699 
5*  .46  9 

5  3.  54  0 
50,524 
50.433 
40.423 
50.497 
51.552 
50, 425 
50.516 
5r  •  5  1  4 
5f. 447 
50 .453 
50.345 
50. 253 
53.537 

50.513 


PH  A  2 

PMfl 

CRLC 

14t.f57 

§P.7§5 

-0,  215 

143.531 

1  13.574 

0.017 

151. 160 

84.737 

0.47? 

154.564 

70.043 

0.443 

150.0?? 

76. 3 CO 

0.263 

149. 106 

90.296 

0.98? 

151.363 

84  ,?51 

-C .349 

154. 1 76 

91 . 615 

0.27ft 

144.969 

7?. 433 

0.  872 

149.739 

83  .618 

0,0  17 

13?. 636 

73.619 

0.  |66 

139.417 

96.? 35 

1.291 

131. 638 

86. 874 

0.027 

153. 539 

77.  3C  3 

-0.49? 

I’?.  40? 

94 74 

1.978 

1  24.  78?. 

90.1  3  5 

**  .  9Q' 

1 ?  0 , 96  8 

97. r 1 4 

7  ,7 

116.547 

193.419 

0 , 9  A*. 

1 2 "  . 5?  8 

8T.  51  5 

-'•.4  9? 

109 . 86 ) 

104.710 

1.141 

1  0  9  .  p  1  > 

1)?.7?2 

1  ,7AT 

92.064 

41.  1  *’? 

(.60? 

1°1  .C  28 

IO6.1 75 

n.71<. 

79.  ]4  ■» 

97.471 

0 .70J 

73.061 

96.?3ft 

0.4?P 

69.939 

1 18.268 

0 . 4  ?  7 

PI. 464 

1 )8.?66 

1.131 

60.7?  7 

91.110 

r.750 

66,Oq> 

1  n5 . 8  1  1 

1.471 

f  2  .  1  1 

94. 041 

0. 888 

48.5  1  3 

1  18.856 

1,356 

61.619 

116.30  3 

7  .597 

64,r*3 

84. ?4? 

0,41  3 

56 

97  .03  1 

1  .40  i 

41  .946 

f'Q.O  7') 

1 . 5  ?  7 

CORRECTED  R,  EL  AND  AZ  FOR  INPUT  TO 
NRTPOD  {ALSO  OUTPUT  ON  PUNCHED  CARDS) 
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APPENDIX  C 


SUBROUTINE  OAKOS  FROGRAM  LISTING 


u  1 


20 


60 


61 


bU 


62 

63 


70 


72 


DIMENSION  DFPG (2)  ,  SUMS  X  ( 4  30)  , I  DENT { 1 5) , I HM ( 2)  ,UPLOG(4,30)  ,  0**00010 

1  SUM  A  (30)  , SUME  (30)  , PIT AG (4)  , PI OT AG  (4)  , MS  EC ( 1 20 0)  ,  IS  EC (  1 20 0)  OAK  000 20 

COMMCN/OAK/TIM  (1200)  , TEN (1200)  ,CRLC{1200)  ,GTNAX(120u)  ,AVGSX  ( 4 , 30)  , 0 A K COO  3 0 

1  DEL  (30)  ,  R  R ANG (12  00)  , R  F  LE V  (1200)  ,AZCR(  1200)  ,DA7{30)  ,IHRX(1200)  ,  OAK00040 

2  IK  I N  X  (120 C)  OAKOOO.SO 

COMMCN/BDCOMT/TIKES (300)  , XSPHA (4,30,300) , R A  NGK  N  (300)  ,ALSAV(  300)  ,  OAKCC060 

1A7I  (300) ,  ELE (300)  ,IRGA(30)  ,IPOL(4)  ,NPOL  OAK00070 


TOUBIE  PRECISION  AVGTfl ,SEC,SUNTH,T1 ,T2, TINES, TOTIM,TSTART,T. STOP,  OAKOOORO 
1TSV,ZSEC,ZSEC1,ZSFC2,Tir,TLIFT,TEM,TM2,PRAC 
TOTIM(IH,IM,SFC)  = C FIC AT ( 6C * ( 60* I  fU  IM) )  *SEC 
DATA  IM  IL/1000000/ 

DATA  10  RBTI/O/ 

DATA  F A  DA  R/' ALT ' / 

DATA  SUME/30*0.0/ 

DATA  SUMA/30*0. 0/ 

DATA  SU MS X/ 120*0.0/ 

DATA  UP  LOG/ 120*0.0/ 

DATA  FLOTAG/' UNCORPECTFD  V-IC'/ 

DATA  FLTAG/'COFHFCTEC  V-LC  •/ 

CALL  REREAD  (90,630) 

TSTO  P  =  0 . 0 
IAG A IN=  0 
SIIMTM=0 . 0 
ICAR  C  =  0 
N  PT  S  =  0 
I ND  E  X  =  0 


OAK00090 
OAKOO  100 
OAKOO 1  10 
OAK  0 0 1 2  0 
0  A  *'  0  0  1  3  0 
CAK001U0 
0  A  K  0  0  1 S  0 
OA  K  00 1 b  0 
OAKOO 170 
OAK  OO 1 30 
OA* 00  190 
CAKO  0200 
OA  K  002  1 0 
OA»’002  20 
0  A  K  0  0  2  3  0 
OAK  0024  0 
OAKC0250 
OA»00260 


READ  (5,20, END=90 1) I  DENT 
FORMAT ( 1 5A4) 

RE  AD (5,  6  0) IM 1 , IM 1 ,ZSEC1,IH2,IN2 , ZS EC2 , NPG, I  NT  APS , I  PAT 
1 ,TSG,TAVG,TSKIF,IC ARD,  INRT 
FORMAT (2  (213, F7. 3)  ,4X,4I5,2F10. 3,215) 

IF  (INRT  .GT.O)  READ  (5,61)  HEAR  ,  I  MONTH,  I  DAY 
FORMAT  ( 3  IS) 

IPOL  (1)  =1 
IPO L  (2)  -=3 
UCL  (3)  =4 
TPOL  (4)  =0 


OAKOO 2 70 
0 A  K  00  28  0 
OAK002°0 
OAKOO 390 
0  A  K  0  0  3  1  0 
OAKCO320 
OAKC 0 330 
OAKOO  340 
OAKOO  ISO 
OAKCO  3  b  0 
0  A  k  0  J  3  7  0 


ISTGT  =  1 


0  A  K  0  0  3  8  0 


WFITE (6, 64) 

FORMAT (SX, 'THESE  ARE  YCUR  INFUT  CARDS') 

WFITF (6,6  2) IH1,IM1,ZSEC1,IH2,IM2,ZSEC2,NRG, INTARG, IP  AT, IPOL, 
1ISG,TAVG,TSKIP 

FOP  MAT  (2  (21  3,  F7. 1)  ,4X,fll3,2F10.  3) 

WFITE  (6,63)  ISTGT 
FORMAT (15) 

CALI  STCIDV (IDE NT,  55,0) 

IF  (NRG. FQ.O) GO  TO  550 

IF  (INTARG. EQ.O) GO  TO  560 

IF(  (IPAT.  EC.O)  .OR.  (IPAT.GT. 3) ) GO  TO  570 

IF  (TAVG.EQ.0.0)  GO  TO  590 

TSTART=TOTIM (IH1,IM1,ZSEC1) 

IF( (TSTAPT.GT.TSTOP) .AND. (IAGAIN.nE.44))GO  TO  72 

CALL  REW 

IAGAIN=1 

I  STOP  -TOTlM  (IH2,IM2,7.SEC2) 


OAKCO  3  )0 
OAK004U0 
OAK004 1 0 
OAKOO 420 
OAFOOU  30 

OAKP0440 
OAK  0045  0 
OAK00460 
OAKOO 470 
CAKOOUSO 
OA  K  004  9  0 
OAK  00500 
OAK005 10 
OAKCOSn 
CAKC05 JO 
OAK00S40 
OAfC05r0 
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CCU  NT -0 
INT=  0 
T1=T STA  RT 
1 2=T  1 tT  A VG 

100  CALL  ALBEAD(TSTART, TSTO P.TLIPT, I  NT ARG , IPAT, 2,  NPTS,NFPG, I AG AIN, 
INRG.ISG) 

IF  (IAGAIN  .FQ.55) GO  TC  901 
IF(IAGAIN.EQ.44) GC  TO  5 10 
1 10  IF (NPTS.EQ.Q)GO  TC  91 

IF  (INT.EQ.O) WRITE (6,1 40) IDE  AT, I NTA  RG 
140  PORN  AT  { '  1  '  ,  10X '  ALTAIR  CAK  VERSION  27  JAN  1971  '//31X, 

1  15A4/31X, 'TARGET  NUNBER  =  ',15//) 

IF ( (INT.EQ.O)  .AND.  (ICARD.GT.O) ) WRITE {7, 150)  IDENT 
1C0  FCRfl  AT ( 1 5  A4) 

INT=  1 

DO  220  i =  1 , NPT  S 

160  IF (TINES  (I)  .GT. T2) GO  TO  240 
IF (TI.GT.TIHES(I)  )GO  TO  220 
i  m  t = i 

DO  210  N=  1 , NPOL 
DO  200  K= 1 , NRG 
EXTEN=(XSPHA(N,K,I)/10.) 

IF (EXTEN.GT.75.0) GO  TO  245 

XSPHA  (N,K,I)=10.**EXTEN 

SDflSX  (N,K) =SUHSX (N ,K)  ♦XSPHA (N, A, I) 

IF (N.GT.  1 ) GO  TO  2C0 

1H5  XSPHC=XSPHA(2,K,It150)  -XSPH  A  ( 1 ,  R  ,  I  ♦  15  0) 

IF  (AES  (XSPHQ) . LT. 180.) GO  TO  190 

IF(XSPHQ.GT.O.O) X SPHA ( 1 , T , I ♦ 1 50)  =X SPH A { 1 , K, 1+ 150)  +360 . 

IF  (XSPHQ.  LT.0.0)  XSPH  A  ( 2 ,  K  ,  I  ♦  1  50)  =XSPHA  (2, K, It  150)  060. 

GO  TO  185 

190  XSPHQ  =  ABS  (XSPHQ) 

SUHA  ( K)  =SUHA  ( K)  ♦XSTPO 

19  1  XSPHX=XSPHA (3, K, It  150) -XSPHA ( ' , K , I ♦  150) 

IF  (ABS  (XSPHX)  .L7.  IPO.)  GC  TO  '(92 

IF(XSPHX.GT.O.O)  XSPPA  ( 1 ,  *  ,I'-150)  =XS  PH  A  { 1 ,  K ,  It  150)  060. 

IF (XSPHX. LT.0.0)  XSPHA (3, A ,  Jt150)  =X SPHA  (3,K, It  150) ♦ 360. 

GO  TO  191 

192  XSPHX  ^ABS  (XSPHX) 

SUI1E  (K)  =  SUHE  (K)  tXSPHX 
200  CONTINUE 
210  CONTINUE 

CCU  NT  =CC UNTt 1 
220  CONTINUE 

IF(IAGAIN.EQ.O)GC  TC  240 

NPTS=0 

GO  TO  100 

240  IF (COUNT. NE.O. 0) GO  TO  280 
245  WRITE (6,260)T2,EXTFN 

260  FOFN  AT (/25X • AT  TIHE  =  • P 1 2. 4 , 2X • TH ERE  IS  A  TINE  GAP  CP  BAD  DATA 
1  EXTEN  =  • ,F10.4) 

GO  TO  440 

280  DO  150  H=1,NPCL 
DO  340  J=  1 , NRG 
I F ( H . GT . 1 ) GO  TO  290 


OAK0056  0 
OAKC0570 
CAKC0580 

CAK00590 

CAK00600 

OAK00610 

CAKC0620 

OAK00630 

0  A  006 4  0 

OAKC0650 

OAK00660 

OAKC0670 

CAKC06P0 

OAKC0690 

OAK09700 

OAKC0710 

OAK C 0720 

OAK  097 10 

OAKC0740 

CAKC0750 

OAKC0760 

OAKCD770 

OAK00780 

OAKC0790 

CAKC0800 

CAKC08 10 

OAK00920 

C AKCO0  30 

OAK00P40 

OAKC0R50 

0AKC0860 

OAK00870 

OAKCC>80 

0  A  E  0  C 1 . 9  0 

OAK  C9900 

0AKC9910 

OA  K  CO 9  20 

OAKC0930 

OAK  C094C 

OAKC0950 

OAKC0960 

0  A  K  CO  97  0 

OAKf 0980 

OAK  0099  0 

OAKOIOOO 

0  A  K  0  10  10 

C  AK01 J20 

0  A  K  0 1  0  3  0 

r  A  K  o  1  0  4  0 

0 AK01050 

OAK01060 

OAKC1070 

OAK01C80 

CAK 0  1Q9Q 

OAKOI 100 
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daz(»:;=suba(J)  /count 

10„  DEL (J)=SUHE  (J) /COUNT 

*  J^S?T(n^,=S0"SX(,,'J,/CCUNT 

IF  (UILOG  (11,. 1)  .ST. 0.0)  GO  TO  100 
AVGSX(H,.1)=9q.q  U  ,U0 

GO  TO  320 

SUNA  (J)  =0.0 
SUNE  (J)  =0.0 
1  CONTINUE 
ISO  CONTINUE 

NIDP7= {  (COUNT* 1)  / 2 J  ♦!  NPT 
AVO.  EL  =  EI  E  (HIDPT) 

A VGAZ*AZI  (HIDPT) 

AVGTH=TIHES(HlnPT) 

AVGRG=R  ANGKH (HIDPT) 

A VGA1=ALSAV  (HIDPT) 
kcunt=ccunt 

^  I  F  ( I  STO.T.  EC.  0)  GO  TO  440 

1ICARC,TLIFT)Vr,AT,'r^TGT'AVRAZ'AVGEL'rGAT»AVi;T'’,  AVGRG.NRG,  INDEX, 
CO  II  NT  =  0 
T  1=T2*TSKIP 
T2=TUTAVG 
SUN  T 1  =  0 . 0 
SUNRO,  =  0. 0 
SUNAL=0.0 

JIlJJ-iM””*””**60  To  160 

I  F  ( T  2.  G  T  .  TSTOP)  GO  TO  *10 
DO  460  K  =  I  NET, NPTS 
KNPT=K 

w»n;ii'TI"ES"”  G0  T°  ',s° 
ne=n  fts-knpt* 1 

DC  501  N  =  1 ,  [ID 

KL=KNPT*N-1 

TINES  (N)  =  TI HES  ( NL) 

RANGHn (n) =RANGKN (NL) 

ALS A  V (K) = ALSAV (NL) 

AZI  (N)  =  A Z  I  (NL) 

EIE(N)  =  ELE  (NL) 

DO  490  K= 1 , NPC  L 

DO  SCO  L=  1, NRG 

XSPH A (K, L,N) = IS PH A ( E,  L , NL) 

CONTINUE  '  '  ' 

CONTINUE 
CONTINUE 
NET  S  =  ND 

IF(I AGAIN. NE.O) GO  TO  100 

s:s;«!;;as,„.ss . . . 


OAKOI 1 10 
OAK01 120 
OAKOI  no 
OAKC  1  140 
OAKOI ISO 
OAKOI  ISO 
OAKC 1 170 
OAKOI ISO 
OAKOI  no 
OAKC1200 
CAK01210 
OAKO  1220 
OAE31230 
OAK01240 
OAKO 1 2  SO 
OAKO  1260 
C A  K  0  1  27  0 
OAK01230 
OAK01290 
OAKOI 300 
oakoi  no 
OAKO  1320 
0  A  K  0  1  n  0 
OAKC  1  340 
OAKOI  )50 
OAKOI  360 
OAKOI  370 
OAKOI  no 
OAKOI 390 
OAKC  140  0 
OAK01410 
OAKO 1420 
OAKO 1430 
OAKO  1440 
OAK01450 
OAKO  1460 
OAKC1470 
OAK014R0 
OAKO  1490 
CAK 01500 
OAK  0  IS  1  0 
OAKO  1S.?0 
OAKC1S30 
OAK01S40 
OAKC1560 
OAK01560 
OAK  C 1 57  0 
OAKC 1580 
OAK  01590 
0  A  K  C  1  6  0  0 
OAK  016  10 
OAKC  1620 
CAKC  Io30 
OAX01640 
CAKC16S0 
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541  FORMAT (40  X, 3F10. 3) 

IF (INRT.EQ.O) GO  TO  5411 
DO  54100  P1=1#  INDEX 
TM2=IDINT  {TFfl  { M)  ) 

ISEC  (M) =TN2 

FPAC=TEM  (M)  -TM2 

IFR  AC=I  DINT  (FR  AC*I  fl  IL) 

MSEC  (M)  =DPIOAT  (TFR AC)  /  1 0 .  ♦.  5 
54 10  C  CONTINUE 


OAK01660 
OAK01670 
OAK01680 
OAK  0 1 59  0 
OAKC1700 
OAK01710 
OAK  0 1 7 2 0 
OAK01730 


lisEcm f  J H *T •  IH» * C X1 , iminx  (i)  ,isec (i) ,  o!kci75o 

541C  FORMAT^ ,  Iro  V  l  U  {l)  '  (I>  ,1=1,  INDEX)  OAK01760 

5.,,  OJJ01770 

HPITE  (6,542) INDEX  OAKC17B0 

5<J"'  OORTOT900//I5''CORRFC7ION  CARDS  HAVE  3EEN  PUNCHED’>  OAK01800 

550  HPITE (6 ,555)  OAKC1810 

SC,S  GORTOT9C02X'A  ZERC  VAl"E  “AS  INP'JT  F°R  NRG  TRIS  IS  *  N0  N0’>  OAK01310 

560  HPITE  (6 ,565)  OAKC1840 

S6510r3Is(GOO?iJCZARSlJA)tUE  “AS  INP°T  P0P  T"E  TARGET  THE  CAT*  EDI-SI?oJUS 

GO  TC  900  OAK  C 1 870 

570  HPITE (6 ,575) IPAT  OAK01880 

™  1^°?  MAS  INP0T  F0R  IPAT  ™E  ™LY  LEGAL  VAOAKCIIOO 


1UIES  FOR  IPAT  ARE  1,2,3')  L 

GO  TO  900 
590  HPITE (6 ,595) 

<’<,rj  in°«r  AT(//2X'A  ZER0  VALUE  F0B  TINC  CAN  NOT  WORK  IT  HILL  BE  SET  TO 
10. 0>  SECONDS  AND  THE  PROG  HA  F  HILL  CONTINUE') 

TAVG=0.C5 
GO  TO  70 
900  I AG  AIN=  9  9 
GO  TO  4  1 
901  CALL  PLTND 
RETUFN 
END 


OAE01910 

OAK01920 
OAKC1910 
OAK01940 
CAK01950 
OAK  C 1 96  0 
OAKC  1970 
OAKC1980 
OAK01990 
O  A  K  0  20  0  0 
OAK020 10 
OA  K  0202  0 
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APPENDIX  D 


*J  \ 


SUBROUTINE  OAKOS  FLOW  DIAGRAM 
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\MRITE  INVALID/ 
VNPJTDATAf 
\  ,v\S6  / 


APPENDIX  E 

SUBROUTINE  GLMP  PROGRAM  LISTING 


S OB  BOUT IN E  GUI  P ( AVG AL, ISTCT , AVG A Z, AVGFL, IGAT,  AVGTN, AVGRG , NRG 
1» INDEX, ICAHD.TLIFT) 

DIHENSION  IH1  (2) 

DOUBLE  PRECISION  AVGT H, TINES, ZSEC, TIN, TEN, TLIET 

COMNCN/OAK/TIH  (1200)  .TEN (1200)  ,CRLC  (1200)  ,GTNAX  ( I  2 0 0 )  ,  AVGSX  (4,  30) 

1  DEL (30)  ,RHANG (1200)  ,FELEV  (1200)  ,AZCR  (1200) ,DAZ (30) ,IHRX  (1200) , 

2  INI  NX (1200) 

CON  FiCN/RDCO  NT/TINES  (30  0)  ,  XS  PH  A  (  4 , 30 , 3  00)  ,  R  A  NGKN  (  300)  ,  A  LS  A  V  ( 300)  , 
1AZI  (300)  ,ELE  (300)  ,  IBG  A  (30)  ,IPOL(U)  ,  NPOL 
EQUIVALENCE  (I  Hfl  (1 )  ,IHR)  ,  (I BN  (2)  ,ININ) 

G  K A  X  =  - 100  0.0 
DO  20  I=ISTGT, NRG 
IF(AVGSX  (1,1)  .LE.GMAX)  GO  TO  20 
GNAX  =  AVGSX  (1,1) 

IGAT  =  IRGA  (I) 

NN=1 

20  CONTINUE 

CALL  TSPLIT (AVGTN, IHN, ZSEC) 

PF1A7  =  DAZ  (NN) 

PHEl=DEL (NN) 

AZ  =  A  VGS  X  ( 2 , N N) 

EL  =  A  VGS  X ( 3,NN) 

CALL  SUBOAK  (IHR, ININ, ZS EC, GNAX,AVGAL, IG AT, A VGTN , A VGAZ , A VG EL, PHAZ, 
1HFL, AZ, EL, INDEX, AVGRG, ICARC ,TLIPT) 

RFTUFN 

END 


PEAC0010 
DE A  C  00S  0 


PE  A  00  Of  0 
PEA  0007  0 
PFAC00P0 
PE  A  0009  0 
PEA00100 
PEA031  10 
PEAOO 120 
PEA001 30 
PE  A  00  1  U  0 
PE  A  00 1 S  0 
PFAOOIf  0 
PFAOO  170 
PFAOO  IPO 
PE  A  00  1  9  0 
PFA00200 
PPFA002 1 0 

PEA  00.3  3  0 
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GMAX  =  CROSS 
SECTION)  (I) 
GATE  NJO.  = 

1  IRGA(X)  I 


SAVE 
AZ  fDGEM) 
E  U  (DBSNf 


'suboakX 

COMPUTE 
OFF-AXIS 
v  CORRECT  I0WS 


RETURN 


J 
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APPENDIX  G 

SUBROUTINE  TSPLIT  PROGRAM  LISTING 


SUBROUTINE  TSPLIT  (A  VGTB  ,1 H B  ,TRU  N) 

TSPC001 0 

DIBENSICN  IHB  (2)  ,  DIVIDE  (2) 

TSP00020 

DOUBLE  PRECISION  AVGTfl  ,TRUN 

TSPC0030 

DATA  DIVIDE/3600. ,6C./ 

TSP00040 

TRUN= AVGTH 

TSPC00S0 

DO  20  1=1,2 

TSP10060 

I H PI  (I)  =TRUN/DIVIPE(I) 

TSP00070 

TRUN=TR UN-FLO  AT {IHB (I)  ) ♦ C I V  IDE  { I) 

TSPC00P0 

CONTINUE 

TSPOOUDO 

PETUEN 

T  S  P  0  0  1  0  0 

END 

TSPOO 1 1 0 

40 


APPENDIX  H 

SUBROUTINE  SUBOAK  PROGRAM  LISTING 


SUBROUTINE  SUBOAK (THB, IMIN, ZSEC, GM A  X , A VG AL , I G AT, A VGTM , AZ E , ELKX , 
1PHAZ,PHEL, AZ, EL, INDEX, AVGRG ,TCARD, TLIFT) 

DIMENSION  A  ( 2  4 )  ,  B  (24)  ,C(27)  ,  D  (27 )  ,  Q  (27)  ,R(16)  ,S(16) 

DOUBLE  PRECISION  ZS EC , TIM , AVGTN , TEM , TL IFT 

COM MON/O A K/TIM  (1200) ,TEM ( 1 200)  ,C BLC  ( 1 200) , GTM AX ( 1 200)  f AVGSX  (4, 30) 

1  DEL (3  0)  , RB ANG  (  1200) , RELEV (1200)  , AZCR  ( 1  200) , DAZ ( 30) , I  HEX (  1 200)  , 
2IMINX  (1  200) 

DATA  DR, NIT, NN IT, NITT/.O  174533,  24, 27,  27/ 

DATA  XKMFT/. 0003040/ 

DATA  A/- 20. 6, -27. 3, -23.  2, -20. 1  ,-17. 0,-12. 8, -10.  4, -0.  4  , -6.  2, -4. 6, 
1-3.  2, -1.7, -0.7, 1.1, 1.7, 2. 9, 4. 1,5. 1,7. 0,8. 5, 7. 5, 10. 3, 10. 8, 12. 4/ 
DATA  B/C.  0,0. 03,  0.  1, 0.17  , 0.  24, 0.  37, 0.59, 0.75, 0.90, 1.05, 1.26, 1  .  34, 

1  1. 48, 1.61  ,1.75,  1.08, 2. 01, 2.  14,  2.  28, 2.  39, 2. 61,  2.  74, 2.  8  1,2. 9 3/ 

DATA  C/- 34. 1,-29. 3, -24. 2, -19. 9, -17. 4, -13. 8, -11. 3, -8. 9, -6. 8, -5.0, 
1-4. 0,-2. 3, -1.1, -0.1, 1.4, 2. 6, 3. 5, 4. 8, 5. 1,7. 9, 9. 6, 11. 3, 12. 7, 14. 2, 
217.  2,18.4,19.7/ 

DATA  D/C. 0,0. 07, 0.1 4, 0.21, 0.28, 0.42, 0.5b, 0.71, 0.05, 1.0, 1.14, 1.29, 

1  1. 44, 1.57, 1.71, 1.85, 1.98,  2.  12, 2.  25,  2.  39, 2. 53,  2.  66, 2. 79, 2. 92, 3. 04, 
23. 17,3. 28/ 

DATA  0/0.0,0.0,0.0,c.0,0.1,0.1,0.4,1.3,2.2,2.9,3.6,5.1,6.3,7.5, 

19.  1  ,  1  0.  9,  12. 5, 15.  0,  17. 5, 1  9.  6, 24.  0, 27.  9,  31. 7, 36.  9,  43.  0,48. 7, 52.  3/ 
DATA  f)/0.  0,0.  02  5, 0.10, 0.17, 0.24,  0.2  9,  0.  35, 0.50, 0.61, 0.65, 0.75, 

1C. 84,0.56,1.00, 1.07,1.  14/ 

DATA  S/ 0.0, 0.25,  1.00,  1.25, 2. 00, 3. 00, 3. 60,  5.  25, 7. 00, 8.  10,  11. 50,  14.' 
1, 18.40,23.70,30.25,40. 25/ 

INuEX=INDEX+1 

DEL  AZ- A7.-GMAX 

DFLEL  =  EI.-GMAX 

XD=99.9 

DO  20  J=  1 , N IT 

IF (DFLAZ.GT.A  (J) ) GO  TO  20 

IF (DELAZ. EQ. A  (J) ) GO  TO  15 

IF ( J. EO. 1) GO  TO  14 

XA  =  A  (J)  -A  (J-1) 

XE=  B  ( J)  -  E  (J-1) 

7C=DELAZ-A  (J-1) 

XD=  (  (XC/XA)  *XB)  +  B  (J-1) 

GO  TO  25 

14  XD=  0. 0 
GO  TO  25 

15  XD  =  B  ( J) 

GO  TO  25 

20  CONTINUE 
25  AF=DR*PH AZ 
AD=CCS ( AB) 

IF (AD.LT. 0. C) XD=-XD 
XH=99.9 

DO  30  K= 1 , NNIT 
I F ( DELEL. GT. C  (K) ) GO  TO  30 
IF(DELEL.EQ.C(K) ) GO  TO  32 
I F ( K . EQ . 1 ) GO  TO  31 
X E=C  (K)  -C  (K-1) 

XF=  D  (K)  -D  (K-1) 

XG  =  DEIEL-C  (K-1) 

XH=  (  (XG/XE)  *XF)  *D  (K-1) 


SUB090  1  0 


SUB00100 
SMRCC 1 20 
SUB 001 30 
SUB00140 

SO  900160 
SOB  00  17  0 
SOP  001 80 
S  U  B  0  0  1 9  0 
SU90020C 

59000220 

SUP00230 

SUBC0240 

SOB  002 5  0 
SUB00260 


SUB  0029  0 
S  U  B  C  0  3  1 0 
SUBC0320 

SUB  00  3 4 0 
SUBOO  3C0 
SUB00360 

SUBOO 7° 0 


no  TO  35 
31  X  H=  0 . C 
GO  TC  35 
3  2  XH=D(J> 

3C  TO  35 
30  CONTINUE 
35  AF=DR*pnEL 
AF-CCS ( A  E ) 

IF  (AF. LT. 0. 0) X H=- X H 
CFL  EX=COS  (FLEX) 

DFLDG=SCPT( (XH**2)  ♦  ( (  (CEL E  X )  *  (X  D)  )  * *2)  ) 

X 0=0 9.0 

DC  4 C  L  =  1  ,  N  ITT 

IF ( DELDG . GT. D (L)  )  GO  TO  40 

IF (DELDG. EQ.D  (L) ) GO  TO  45 

X0=  (  (DELDG- C  (L-1)  )/(C  (L)-D  (L-1))  )  •  ( Q  (  L)  -Q  (L- 1 )  )  (L-  1 ) 

GO  TC  4 1 
4C  XJ-Q  (L) 

GO  TC  41 

40  CONTINUE 

41  X  UHF= c0.0 
NTX=  16 

DO  5  C  1=1,  NIX 

IF (DELDG. GT. R  (1) ) GO  TO  50 

IF (DELDG. EO. P  (1)  ) GO  TO  51 

XGHF=  (  (DELDG-R  (fl-  1 )  )  /  (P  (R)  -  P  (R- 1 )  )  )  »  (S  (R)  -S  (P-1)  )  »S  (P-1) 

GC  TC  52 
c  1  XFHF=S(R) 

GO  TO  5? 
c  0  CONTINUE 
*2  AT  E=  A  2S/D  P 
ELE  X  =  EL  EX/PP 
I  HP  X  (INDEX)  =IHP 
IPXNX (INDEX)  =  T I N 
TIE.  (INPFX) =AVGTR-TLIFT 
TFE  (INDEX)=ZSJC 
GTRAX (INDEX) =  GRAX 
CFLC (INDEX)  =GF  A  X ♦ XD 
AVGPG=A VGPG*  (  ( IGAT»30. 0)  / 1000.  0) 

KF= AVGPG/XKNFT 
E=EL  EX 

CALL  FEFC  (E , PS , DF E  ,  CPF ) 

PFANG  (INDEX)  =AVGPG-  (DRR*XK*,FT) 

PEL S V (INDEX) =ELEX-DEE*XH 
AZCP (INDEX) = AZE»XC 
XF=XDHF 

IF (ICAR D. EC. 0) GO  TO  54 
•FITE (7,53) A  VGTH , XK, X J 
c’  Ff 3"AT (3F 10. 3) 

54  T. F  (INDEX.  GT.  1)  GO  TO  70 
»'FITE  (6,55) 

5 c  fOPRA"r(52X' DELTA  DELTA  DELTA  DELTA  DELTA') 

«' FITE  (6,60) 

61  FOPNAT  (SX'TIFE  AZ  FI  V-IC(EB)  AZ  (DB),  F.L(D3)  AZ  DEG  EL 

1  DEGHEES  VHP  (CB)  OBP(DP)  PHAZ  PR  EL  CPLC') 

10  y FITE (6, 7 5) IHR,IHIN, ZSEC,AZE, ELEI,G RAX, AZ, EL, XD,IH, DELDG, XJ,XN 
1Z,  PFEl.CRLC  (INDEX) 

7  5  FOPNAT (2I3,F7. 3,2F7. 2,F6. 2, 3XF6.  2, 1XF6.2,2XF6. 2, 1XF6.  2,3 (2XF6. 
1 :  F  1 C  .  3 ,  F  1  0.3) 

S  FT  U  F  N 
F  N  2 


SDPC0390 
50500400 
GUB004  10 
SUB00420 

SUP  204  30 
SUP 004  4  0 


SUBCC4P0 
SUB004O U 
SUB  C  0  5  1  0 
SUB  CO  5 1 0 
SUB  00530 

SUP00550 
S  U  P  0  0  5  6  0 

SUB006  /O 


.SUB0J6OL’ 


SUPC0620 


SUB  C  '64  0 
SU5C065C 


SUBC0700 
SOB'  ,710 

DEGSUBC0730 


PHA 
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APPENDIX  J 

SUBROUTINE  SUBOAK  FLOW  DIAGRAM 


% 


THE  TABLES  CONTAIN: 

A  AAZfDB) 

B  AA2 
C  AELfDB) 

D  A  EL  $  A&" 

Q  A  VHP  RCE(DB) 

R  A& 

S  AllHF  RCb  (D3) 


INDEX=lWDEXtl 

COMPUTES 

AAz^b)+Ael(db; 
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REFC 

GET 

DRR<  DEE 


APPENDIX  K 

SUBROUTINE  REFC  PROGRAM  LISTING 

SUBROUTINE  REFC(F,F,CEE,CRR)  VERSION  6/16/70 

DIMENSION  DE  (16,8)  ,CR  (16,8)  ,ED<16)  ,RD(8) 

CATA  CE/O.C  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  , 

10.0  ,0.0  ,0.0  ,0.0  ,C.O  ,0.0  ,0.0  ,0.0  ,0.0313, 

20.030  3,0.  0292,0.0  28  7,0.028  2,0.  02  72,0.  0262, 0. 0 25 3 , 0. 024  3 , 0 . 022 3 , 
3C. 0214,0. Cl  95, 0.0171, 0.0135,0.  00 75,0. 0  ,0.0937,0.0848,0.0770, 

40. 0732, C.C6;4,0.C627,0.C57 1,0. 0522,0. 0480, 0.0412, 0.0385,0. 03 37, 
50.0  27  8,0.0205,0.0105,0.0  ,0.  18 50,0.  1520, 0.1 250, 0.1  1  40, 0.  1050, 

60. 0 9C 4, 0.C795.0.C708, 0.0636, 0.0523,0. 04 78, 0.0405, 0.0323, 0.02 29, 
70. 01  14, 0.  0  ,0.53 10, 0.3070,0.  21 20,0.  19 30, 0.  1 60 0,0. 1280, 0.  1 06  0, 

8  0.0899,0.  0780,0.061  2,0.0550,0.  0455,0.  0354 , 0 . 0246 , 0.  0 1  20, 0 . 0  , 

9  0.7550,0.  3720,0.  2400,0.  20  20,0.  1750, 0.  1370,0. 1  120,0.0  942,0. 091  1  , 

A0. 0631, 0. C566, 0. 0466, C. 0361, 0,02 50, 0.0122, 0.0  ,0.9 120, 0.41 10, 

B0.  2  560,0.  214  0,0.  18  40,0.  14  20,0.  1150,0.  0967,0.0  83  0,0.06  4  3,0.0  575, 
CO.  0472,0.  0365, 0.0 25 2,0.  0122, 0.  0  , 0. 97 00, 0. 4 200  , C.  26  CO, 0.2200, 
DO.  1900,0.  1460,0.  1  17C,0.CS80,0.  0840,0.  0653, 0.0584, 0.0478, 0.0369, 
EO. 0254,0. 0123,0.0  / 


PEFCOO  10 
REFC0020 
PFFC0030 
REFC0040 
REF000  50 
REFC 0060 
PEF00090 
PFP00080 
PEF00090 
HEPOOinO 
PEF001 10 
RFF00120 
REF001 30 
PFF00140 
PFFC015 0 
REFC  1160 
PEF00170 


DATA  CR/  0. 

0,  0. 

0,  0. 

0,  0. 

0,  0.0,  0. 

0,  0. 

0,  0. 

0,  0. 

0,  PEFC01B0 

1  0.0,  0.0, 

0.0, 

0.0, 

0.0, 

0.0,  0.0, 

22.6, 

21.5, 

20.4, 

19.9.PFF00180 

2  19.4,  18.5, 

17.6, 

16.8, 

16.1, 

14.8,  14.2, 

13.2, 

12.  0, 

10.4, 

8.6,PEPC0200 

3  7.7,  67.3, 

57.9, 

50.2, 

47.  0, 

44.  1,  39.3, 

35.4  , 

32.  1  , 

29.3, 

24.8.PEFC0210 

4  22.9,  19.7, 

16.3, 

12.7,, 

9.4, 

(.1,132.0, 

98.5, 

77.4  , 

69.7, 

63. 2, 9EE00220 

5  52.9,  44.7, 

38.4, 

33.4, 

26.4, 

23.9,  20.1, 

16.4, 

12.7, 

.4, 

9.1, PFF032  30 

63 40. C, 167.0, 

103.0, 

R6.1, 

73.4, 

56.7,  46.2, 

38.9, 

33.  6, 

26.4, 

24.0, 9EFC0240 

7  20 . 2,  16.4, 

12.8, 

9.5. 

8.2, 

405.0,170.0, 

104.0, 

86.  3, 

73.6, 

56. 3, PEFC0250 

8  46.3,  38.9, 

33.7, 

26. 5, 

24.1, 

20.3,  16.5, 

12.8, 

9-5, 

8.2, 

421.0, PFFC0260 

9171.0,104.0, 

86.6, 

73.9, 

'3.1, 

46.4,  39. C, 

33.8, 

26.8, 

24.  3, 

20.5.PEFC3270 

A  16.6,  13.0, 

9.  8, 

8.4, 

466.0, 

172.0, 105.0, 

97.4, 

74.  0, 

58,0, 

46.6 , PEF00280 

B  19.2,  34.0, 

27.0, 

24.6, 

20.7, 

16.7,  13.0, 

10.0, 

8.  4/ 

9EF00290 

EATA  ED , RTDEG/0. 0  1 

,2.0,4 

.0,5.0 

,6. 0,8.0,  10. 

0,12.0 

,14.0, 

18. ,20 

.,  PEFOiPOO 

1  24.  ,  30.  ,4  0. ,60. ,9  0. ,57.  2657  8/  REF 00  3 10 

DATA  RD/O. Cl, 10. ,30. ,60. ,200. ,400. ,  1000. ,2000./  EFFP032P 

IF(P.LE.O.O)GO  TO  300  REF00330 

90=9/6080.27  PEF00340 

DO  ICC  IED=2, 15  p^FCOIFO 

I  =  1 7 -I E  D  9EFC0360 

IF (E. GE. ED (I) ) 00  TO  120  '’FF07370 

CONTINUE  PPFCC’80 

1-1  PEF.I0  390 

DO  2C0  JRD=2,8  500400 

.i=io-jrd  r-Tcnio 

IF  (RO.GF. RD  { J)  )  GC  TC  220  F9FO0420 

CONTINUE  REF01U30 

J=1  PFFC0440 

IF  (J  .EQ.  8)  GO  TO  340  SSF00450 

7.P=  AlOG  (HG/RD  (J)  )  /AlOG  (90  (J  *1)  /9D(J)  )  9rF01460 

IT  (i. IE. 0.0)  GO  TO  320  P'F0(I470 

ZE=ALOG  (E/ED  (I)  ) /AlCG  (ED  (1  +  1) /ED  (I)  )  i'FrC  ^90 

DF1= ( (DE (1+1 , J)-DE  (I, J) ) *  (1.  -ZR)  ♦  ( DE  ( I,  J+ 1)  -DE  ( I ,  J) )  * Z  9)  *Z  S  RFFC9U90 

DE2  =  (  (DE  (I,  J  +  1)  -DE  (I,  J)  )  *  (1.-ZF)  ♦  (DE  (I*  1,  J*  1)  -DE  (I ,  J  +  1)  )  *7.E)  *7P  PFF0050  0 

DEE=DE1+DE2+DE  (I,  J)  9EFC0C10 

DP1=  (  (DR  (1+1,  Jl-CR  (I,  J)  )  ♦  (1.-ZP)  ♦  (DR  (I,  J+1) -DP  (I,  J)  )  *ZP)  *ZE  9 FF 00620 

DP2=  (  (DR  (I, J+1) -DR  (I,  J)  )  ♦  (1.-ZF)  ♦  (DR  (I+1.J+  1)  -D9  (I , J ♦  1)  )  *7S)  *7. R  R'FOOSID 

0PR=  (DR1+ER2+DR(I,J)  )  D  EFC054  C 

GO  TC  400  PFF0C55C 

DEE-0.  0  ""f"'  ,  R 

DRR=  0. 0  9FFC367  1 

Gi  TC  400  rEFCDCQ0 

DFE=CE(I,J)+(DE(T,J+1)-DE(I,J))*7P  PSFCC5°0 

DRR  =  CR  (1,0)  ♦  (DR  (I,  J  +  1)  -DR  (I,  J)  )  *7R  PEE  006  0  0 

GO  TO  400  RFEC0610 

DFIT*  (E-ED(I) )/(FD (1+ 1 ) -ED (I))  PFFC0620 

DEE=DEIT*  (DE  (T+1,  J)-DE  (I,J)  )  +DE  (I,  J)  FEE  00  6  30 

DRR  =  DEI.T*  (DR  (1+  1 ,  J)  -  DR  (I,  J)  )  +DP  (I ,  J)  PEFC0640 

RET'IFN  9EFC0650 

END  pPFC.1660 
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APPENDIX  L 


SUBROUTINE  ALREAD  PROGRAM  LISTING 

SUBPOUT  IN  If  ALRFAD  (T  START,  TSTOP , TLIFT, I NTA RG,  INPAT,  NOPHA,N?TS,NFP, 
IN EH  FAS, NRG, ISTGAT) 

DIH  ENSION  ALT  (5) , A BT  (3 , 5)  , A  2  (6)  , CALADD (4)  ,CALCON (19) , EL (6) , IHD (13) 
1, IHT  (5)  , INGATE  (5,3)  ,IPRIC5 (5) ,ISENS  (6) ,ISLIDE (5,3) ,ISPAC  (5, 3) , 
2ITARET  (5,3,3)  ,ITARG(5)  ,  IHAVE  (6)  ,  LOC  (5 , 3,  4)  ,  HODE  (5,3)  ,NC(4)  ,R’IX(2)  , 
3NBODES (5,3) ,NPTEST  (2) ,NSABP  (5,3)  ,NSABPT(5)  , POWER (6) , R AD  (6, 5) , 
4RANGE  (5) , VEL(5)  ,IFPAF  (22) 

COflHCN/TREAD/LN,IFLG,IBTRHD,FBTRHD,FBTRHA,FHTRTG,FHTREI, FBTRSP, 
1FNXSBC, F  K ARSG,F8RR1 1, F BGLOT , FBC H AF , FBBSBC, FBAS L P , FBABP  (6)  , FBPHA  (6) 
2, N ABE  (25)  ,  NI  (24)  , IX ( 24)  ,TABP ( 1 28 , 6)  , TPH ( 1 28 , 6)  ,ITEB(2  000) 

cobmcn/tibccb/ibbs,:fs 

COBBON/RDCOBT/TIBES  (300) ,XSPHA  ( 4 , 3 0 , 300) , R A NGKB (300) , ALSAV (300) , 
1A2I  (300)  , ELE(jOO)  , I RG  A  (  30 )  , IFO  L  ( 4)  , NPOL 
EQUIVALENCE  (IHD  ( 1)  , I  DR  EC | ,  (IHD (2)  , LR  EC) ,  (IHD (3)  , ITGTB 1 )  , 

1  (IHD  (4)  ,ITGTH2)  ,  (IHD  (5)  ,IFPG)  ,  (IHD  (6)  ,IPRI)  ,  (IHD  (7)  ,  NELRD)  , 

2  (IHD  ( E )  , BACYBA)  ,  (IHD  (9)  ,  NTARG)  ,  (IHD  (10)  ,NTPSA)  ,  ( I  HD  ( 11 )  ,  N  ’ll  NOR)  , 

3  (IHD  (12)  ,  BICYBA)  ,  (IBD(13)  ,LBICY)  ,  (IF PAR  (1)  ,  FBTRHC)  , 

4  (RAD  (1,  1)  ,  POWER  (1)  )  ,  (RAC  (1,2)  ,AZ(1)  )  ,  (RAD  (1,3)  ,  EL(1)  )  , 

5 (ITARDT  (1,  1,  1)  ,BODB (1, 1) ) , ( ITAR DT (  1 , 1 , 2) ,ISPAC (1, 1))  , 

6 (ITARDT  (1, 1,3)  , N S AB E  ( 1 , 1) ) ,  (IBT(1) ,IGCHG) ,  (IBT (2)  , IGA  IN) , 

7  (IBT  (3)  , NFLS) ,  (IBT (4)  , IS L E A) ,  (IHT (5) , RBP) ,  (ABT (1, 1) , RANGE  ( 1) )  , 

8  (ABT  (1,2)  , VEL  ( 1)  )  ,  (ABT  (1,3)  ,ALT(1)) 

DATA  PCON  ,RKB,VKB/1.0E6,1.873703E-3,4. 4672E-4/ 

DATA  NEX,  NFTEST/0, 1 ,300 ,,150/ 

DOUBLE  PRECISION  DR  A N G , F INT IB , FPG, GBTIBE, PR F , TIB E (6)  , T IBES , TI BOLD , 
1 TLIFT, T STAR T,T STOP, 2TD IF (6) 

IVTEGER*2  iteb 
IF  (NEWPAS.GT. 2) GO  TC  2C00 
IF  (NRG.GT.30) NRG=10 
NPC 1  =  0 
DO  10  1=1,4 

IF  (I POL  (I)  .  NE.O)  NPOL=NFCLM 
10  CCNTINUE 

GBTIBE=TSTART-1.0 
N  PT  S  =  0 
TINC=0. 0 
NTERE=0 

CALL  THEAD  (N EWP A S , S9960) 

DC  60  1=1,6 

IF(IFPAR  (I) .GT.O) GO  TO  60 
WHITE  (6,40)  NABE(I) 

40  FORB  AT  (  '  FORBAT  TABLE  ',A4,'  WAS  NOT  FOUND  -  FUN  ABORTED.') 
NEWPAS=55 
RETURN 

60  CONTINUE 

80  DO  ICC  1=2,19 

CALCCN (I) =GET (FBXSEC,IBTRHD, I) 

100  CONTINUE 
TLIFT=0. 0 

ITL  =  I GET  (FBGLOT, IBTBHD,  1) 

IF  (ITL. NE. 2) GO  TO  1 40 
I RH  S  =  IG  £T  (FBGLOT, IBTPHD, 2) 

I FS  =  I  GET (FBGLOT, I BTRHD, 3) 

CALL  GBTUPK  (TLIFT) 

GO  TO  140 
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120  C  ALL  EREAC(I) 

140  I  HD  (1)  =  IGET  (FHTRHD  ,  IETRHD  ,  1) 

IF  (  (IDREC. LT. 1 27)  . AND.  (IFIG. NE. 3) )  GO  TO  180 
WRITE (6, 160) GHTIHE 

160  FOR  HAT ( '  END  OF  FILE  E  NCCO  NTERED  ON  INPUT  TAPE  AFTER  TIHE  •  F12.4) 
NEWPAS=  44 
RETUFN 

180  IF  (IIREC. NE. 2) GO  TO  120 
TIH01D=GHTIHE 
DO  2C0  1=2,13 

I  HD  (I) =IGET(FHTRHD,IBTRHD,T) 

200  CONTINUE 

I HM  S  =  IC  FT  (FHTRHD, IBTRHD, 14) 

1FS  =  1GET  ( FHTRHD, 3 ETRHD, 15) 

CALL  GHTUPK (GHTIHE) 

IF  (GHTIHE.LT.  (TSTARI-0. 20) ) GO  TO  120 

IF(  (GHTIHE. GT.TI  HO  I,D)  .AND.  ( H  AC  YB  A.  E  C.  0)  .AND.  ( NTDE A . EQ . 0) ) GO  TO  120 
NCO  NT  =  0 

IF  (  (  (CABS  (TIHOLD-GHTIHE)  }  .  GT .  0.  0  1)  .  AND.  (NTDBA.NE.  0)  )  GO  TO  280 
IF  (NTEFR.NE.O)GC  TO  120 
NCC  NT  =  1 
GO  TO  1200 
280  FINTIH=GHTIHE 
TIHOLD  =  GHTIi1E 
N  FP= IFPG+ 1 
DC  300  1=1,5 
N  SA HPT ( I ) =0 
300  CONTINUE 

PFF=FCOK/FLOAT(IPFI) 

TIN  K= 1 . /P  RF 

IF (NBLRD.EC.O)GO  TO  600 

NBLRCC=NELRD 

DO  5C0  1=  1 , NBLRD 

IBTRHA  =  IBTRHD»HACYBA*24*  (I- 1) 

DO  400  Ml,  3 

RAD  (I, K) =GET (FHTRHA.IETRHA, K) 

400  CONTINUE 

IHNS  =  1GET  (FHTRHA,IBTRHA,4) 

I FS  =  IGET  (FHTRHA, IETF  HA, 5) 

CALL  GHTUPK (TIHE (I) ) 

ISENS (I) =IGET(FHTRHA,IBTRHA,6) 

I  HAVE (I)  =  IGET(FHTRHA,IETRHA,7) 

IF  (IHAVE  (I)  . E0.4) I  HAVE  (I) =3 
IF (I  WAVE  (I)  .GT. 3) I  WAVE  (I)  =  0 
500  CONTINUE 

NCH  ANZ= IFPG 

IF (IF PG .EC. 1) NCHANZ=NCHANZ»3 
DO  560  IP=1 ,NPCL 
NC (IF)=NCHANZ*IPCI (IF) 

DO  520  1=1, NOPHA 
I ND  E  X  = (1-1) ♦6*NC(IF)  +  12 
IF(1FFAR (INDEX)  .GT.C) GO  TO  5  20 
WRITE  (6,40)  NAHE  (INDEX) 

NEWPAS=66 

FETUEN 
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5  20  CONTINUE 

KCNCET=3*  (NC  (IP) -  1)  «IHAVE  »1 

CAL  ADC (IP) =C  ALCC  N (KCNGET) -POWER ( 1)  -  FLOAT ( 1 0*1  SEN  5 ( 1 ) ) 

560  CONTINUE 

600  IF (NTARG .EQ.O) GO  TO  120 
ITP I K  =  0 
NBPTST=  0 

DO  9C0  1=1, NTARG 
IBTRTC=IBTRHD+NTEEA*12*<I-1) 

ITARG  (I) =IGET(FHTRTG,IBTRTG,1) 

IPRICF (I) =IGET(FNTRTG,IBTRTG,2) 

NK=  1 

DO  8C0  K=1 ,3 
I K=  E  K  + 1 
N  K  =  I K  ♦  2 
DC  7C0  L=  1  , 3 

ITARDT (I, K,L) = I GE T ( FHT RTG , I ETRT G , I ♦ IK) 

700  CONTINUE 

IF (ISP AC  (I,K) . EQ. 259) ISPAC  (I,K)=-1 
ISP AC  (I , K) =2** ( ISP  AC (I,K)  +NIX (IFPG ♦ 1)  ) 

IF  (ISPAC  (I,K)  .EQ.0)ISPAC(I,K)=1 
NSAHPT  (I)  =  NSAHPT(I)  ♦NSAHE  (I ,  K) 

NHODES(I,K) =0 
DO  760  1=  1,  9 
LOC  ( I , K , L ) =0 
HODUH=NOD  (NODE  (I, K)  ,2) 

IF  (HCDUH. FQ.O) GO  TO  760 
NHODES (I,K) =N NODES (I,K)  *1 
LOC (I, K , l) = NHODES  (I  ,  K) 

760  NODE (I, K) =HODE  (I,K)/2 
780  CONTINUE 

ISIIDE ( I ,  K) =2*NSAHP  (I ,  F )  *NNCDES(I,K) 

NBPTS'l-WtPTST+ISLIDE  (I,K) 

800  CONTINUE 

DO  820  J  =  1 ,  3 

INGATE (I, J)=IGET(FHTRTG,IBTRTG, J*1  1) 

820  CONTINUE 

IF (ITARG  (I)  .EQ. I NTARG)  ITEIF=I 
900  CONTINUE 

NBPTST=NEFTST-HOD  (NBPTST,6) 

NTERR=0 

IF ( ITPIX . GT. 0) GO  TO  960 

WRITE (6 , 920) GHTINE, INT ARG , ( IT A R G (I ) , 1=1, NTARG) 

920  FORNAT ( '  AT  TINE  =  ',F12.9,»  TARGET  ',12,'  IS  NOT  ON  THE  TAPE,  TA 
1RGETS  AVAILABLE  ARF  ’,513) 

IBTRHI=IBTRBD*HICYE A 
NBP= IGET  (FHTRNI ,IBTRNI, 51 
IF (NEP. EQ  .NBPTST) GO  TC  10000 
NTERF= 1 
GO  TO  120 

S 6 0  IF(NHODES(ITPIK, INPAT)  .GE. NECLJGO  TO  1020 
WRITE (6 ,1000) 

16  CO  FOR  HAT ( '  EITHER  PCIABIZATICN  OR  PATTERN  CHOSEN  IS  NOT  AVAILABLE') 
GO  TO  10000 

10  20  IF (ISTGAT.GT.NSAHE (ITE IK , IN  PAT)  )  ISTGAT  =  1 
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IF  (ISTGAT.LT.  1)  ISTGAT=  1 

NST  EST=  NS  ABF  (ITP I F , IN P AT) - 1 STG AT ♦ 1 

NETGAT=ISTGAT 

IPAT=INFAT 

1=1 

10<J0  IRG  A  (I)  =INGATE  (ITPIK, I  PAT)  ♦(NSTGAT-I)  *  ISP  AC  (TTP  IK  ,  IP  AT ) 

NSTG  AT= 1 

1060  IF  (I.GE. NRG) GO  TO  1100 

IF (  ( I ♦ 1 )  . GT . NSTEST) GO  TO  1080 
1  =  1+1 

IPGA  (I) =  I RG  A  (I  -  1)  +ISFAC (ITPIK, IPAT) 

GO  TO  1060 

1080  IFAT=IPAT+1 

IF  ( I P  AT . GT. 3) GO  TO  1100 

IF(NSAHP  (ITPIK, IPAT) . IE. 0) GC  TO  1080 

NSTEST=NSTFST+NSAHP (ITP IF, IPAT) 

IF  ( N  BCD  E  5  (ITFIK,IPAT)  . LT. NPOL) GO  TO  1100 
1  =  1  +  1 

IF ( I.  LF. NRG) GO  TO  1040 
1  =  1-1 

1100  IPAT=INPAT 
N  RG  =  I 

1200  IBTfiHI=IPTPHD+HICYE A 
DO  2200  HIN=1,NHINOE 
IF (BIN. LE.  1) GO  TO  1220 
IBTEBI=IBTRHI+LHICY+NPLS*NEP 

1220  DO  1240  1=1,5 

IBT (I) = IG  El  (FHTRHI,IBTRHI,I) 

1240  CCNTIilUE 

IF (NPP. FQ. NRPTST) GO  TO  1280 
WRITE (6, 1260) GHTlrF,FIN,NBF,NBPTST 

1260  FORBATC  AT  TIBE  =  ',F10.4,',  HINOR  CYCLE  ',11,*  SOHET3ING  IS  SCRE 
1WED  UP  NBP  =  ',15,'  IT  SHOULD  BE  ',14,'  -  SHOOT  H ARTOG ENS IS » ) 

GO  TO  1 3  3  C 

1280  IF(NILS.IE.  1) GO  TO  1340 

WRITE  (6  ,  3"  ’)  NPLS ,  BIN,  N BINGE 

1320  FORBATC  NPLS=  '  ,  1 2,  3  X ,  '  B INO  R  CYCLE  •  ,  1 1,  jX,  1 1  ,  '  HINCRCYCi.ES') 

1330  GBTIBE=GBTIBE+TINC 
GO  TC  120 

1340  TTNC  =  FLC AT  (NPLS) /PPF 
DO  1 5C0  L=1 ,NTARG 
IBTRBT=IETRBI+  (L-1)  *12 
DO  1400  1=1,3 

ABT  (1,1)  =GET  (FHTRBI,IETFFT,I  +  5) 

1 4CC  CONTINUE 

VEL (I)  =  V  E  L (L) *  VK  H 

15C0  CONTINUE 

IBTRBT=I BTRBI+ (ITPIK-  1) *12 
TSD  RE  1= IG  ET  ( FHTRB I ,  IBT R HT  ,  9) 

KSDUB=N  SA  HPT (ITPIK) 

IF  (NRG. LT.NSDUH)  NSDPH=  NRG 
CPANG=RANGE (ITPIK) 

IF (DRANG. GT.0.0) GC  TO  1560 
WRITE (6 , 1540) GHTIB  E , DF  ANG 

1540  FORBATC  AT  TIBE  =  ',F15.4,»  THE  RANGE  =  ',1PE20.6) 


*3S gjpsw 
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DPANG=1 . 0 

I'ifiO  DR  AN  G=  40.  *DLOG10  (DRANG) 

CALIF=DRANG+FLOAT (I GAIN) 

DC  2000  1=1 ,  NPL S 

IF(  (NIN.GT.  1)  .OR.  (I.GT.  1)  )  F  INTI B  =  F INTIN +TINK 

IF (TSTART.GT.FINTIB)GO  TC  2000 

NPTS=NPTS+1 

IF  (NFTS. GT. 1) GO  TO  1600 
1580  TINES  (NPTS)  =FINTIN 
GO  TO  1620 

16CC  IF  (  (I. EQ.  1)  .  AND.  (BIN. EQ. 1) . AND.  (NCONT.EQ.O) ) GO  TO  1580 
TIB  EE  (NPTS)=TINES ( N PTS -  1 ) +T INK 

1620  (ANGKB(NPTS) =RA  NGE (ITPIK) * R K B  +  V  EL ( ITPIK) * ("INES (NPTS) -FINTIB) 
ALSAV  (NFTS) = ALT (ITPIK) 

IF  (NELRD. NE. 0) GO  TO  1630 

INTEFF=NELRCC 

GO  TO  1 6  5  C 

163  0  ZTDIF  (1)  =  CABS (TIBES (NFTS) -TINE  (1) ) 

I  NT  EF  P  =  1 

DO  1640  N AE  =  2 , NBLRD 

ZTDIF (NAE) =C APS (TIBES (NPTS) -TINE  (NAE) ) 

IF (ZTDIF  (NAE) . LE. ZTDIF  (NAE- 1) )  I NTERP=NAF 
1640  CONTINUE 
165C  AZI  (FPTS)  =AZ  (INTER?) 

El F  (NPTS) =EL  (INTER  P) 

ISABFT  =  IBTRHD+ISCFA'»ISCREI  +  (1-1)  *NBP 
NSTEST=NSANP  (ITPIK, INPAT) -ISTGAT+1 
I A  D=  0 

IF  (INPAT. EQ.  1) GO  TO  1680 
JST'-INPAT  -  1 
DO  1660  J=1 ,  JST 
IAC=IAD  +  I SLIDE  (ITP IR , J) 

1660  CONTINUE 
1 6  8  C  L  =  I STGAT- 1 

DO  1800  K= 1 , NRG 
L=L  ♦  1 

IF (K . LE . N STEST) GO  TO  1700 
IAD=I AD ♦ I  SLIDE (ITPIK, I  PAT) 

IPAT=IPAT+1 

NSTEST=NSTEST+NSANP (ITPIK ,1  PAT) 

L=  1 

170C  DO  1780  I  P=  1 , NPOL 

IPIK=ISABPT  +  IAD  +  2* (IOC (ITPIK, I  PAT, I  POL (IP)  )  -1+NHODES (ITPIK, IPAT) 
1*  (L~  1)  ) 

IABF=IGET  (FBTRSP.IPIK,  1) 

IF ( (IANP.LT.1) .OR. (IABP.GT. 128))IANP=1 
XSPHA(IP,K,NPTS)  =TA  BP ( I AB E , NC (I P) ) +  C A  LI B+CAL ADD ( IP) 

1740  IF(NCPHA. EQ. 1) GO  TO  1780 
IPHA=1GET  (FBTRSP.IPIK, 2) 

IF((IPHA.GE.O)  .AND.  (IPHA.LE.127)  )G0  TO  1760 
X5PHA (IP,K,NPTS+150)=0. 0 
GO  TC  1780 

176C  XSPHA  (IE,K,NPTSM50)=TPH(IPHA+1,NC(IP)  ) 

178C  CONTINUE 
1 8C0  CONTINUE 


IP  AT  =  IN  PAT 

IF (TSTOr.LT.TINES(NPTS) ) GO  TO  10000 
IF  (NPTS. LT. NPTEST (NCPHA)  )  GO  TO  2000 
NEHP AS=  99 
RETOFN 

2000  CONTINUE 
2200  CONTINOE 
GC  TO  120 

9960  WRITE  (6 ,9980) 

9980  FORMAT  ( *  THEAD  HAS  DEFAULTED  -  RUN  HAS  BEEN  ABORTED.') 
1C00C  NEWPA£=0 
RETUFN 
iND 


